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Abstract 

We present the multiplicity and pseudorapidity distributions of photons produced 
in Au+Au and Cu+Cu collisions at i/snn = 62.4 and 200 GeV. The photons are 
measured in the region —3.7 < rj < —2.3 using the photon multiplicity detector 
in the STAR experiment at RHIC. The number of photons produced per average 
number of participating nucleon pairs increases with the beam energy and is in- 
dependent of the collision centrality. For collisions with similar average numbers 
of participating nucleons the photon multiplicities arc observed to be similar for 
Au+Au and Cu+Cu collisions at a given beam energy. The ratios of the number 
of charged particles to photons in the measured pseudorapidity range are found to 
be 1.4 ± 0.1 and 1.2 ± 0.1 for ./s^ = 62.4 GeV and 200 GeV, respectively. The 
energy dependence of this ratio could reflect varying contributions from baryons to 
charged particles, while mesons are the dominant contributors to photon production 
in the given kinematic region. The photon pseudorapidity distributions normalized 
by average number of participating nucleon pairs, when plotted as a function of 
f] — Uheam, Sire found to follow a longitudinal scaling independent of centrality and 
colliding ion species at both beam energies. 

Key words: Particle production, photons, forward rapidity, limiting fragmentation 
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1 Introduction 



For high energy heavy-ion coUisions, measurements of particle multiphcity 
provide information on particle production mechanisms [1]. Event-by-event 
fluctuations in the multiplicity of produced particles within a thermodynamic 
picture could be related to matter compressibility [2] . The event-by-event cor- 
relation between photon and charged particle multiplicities can be used to test 
the predictions of formation of disoriented chiral condensates [3]. The varia- 
tion of particle density in pseudorapidity {rj) with collision centrality can shed 
light on the relative contribution of soft and hard (perturbative QCD jets) 
processes in particle production [4]. Multiplicity measurements can provide 
tests of ideas on initial conditions in heavy-ion collisions based on parton sat- 
uration [5] and color glass condensates [6]. Under certain model assumptions, 
the particle density in pseudorapidity can provide information on the initial 
temperature and velocity of sound in the medium [7] . The pseudorapidity dis- 
tributions are found to be sensitive to the effects of re-scattering, hadronic 
final-state interactions, and longitudinal flow [8] . 

Several interesting features of the dependence of particle density in pseudo- 
rapidity have been observed in Au-|-Au collisions from the experiments at 
the Relativistic Heavy- Ion Collider (RHIC). Particle production is found to 
follow a unique, collision energy independent, longitudinal scaling [9] in p+p 
and d+Au, as well as in heavy-ion collisions [10,11]. Such longitudinal scal- 
ing is also found to be independent of collision centrality for photons [12,13]. 
The total charged particle multiplicity (integrated over the full pseudorapidity 
range) per average number of participating nucleon ((A^part)) pair is found to 
be independent of collision centrality [10]. However, at mid-rapidity {\r)\ < 1), 
charged particle multiplicity per (A^part) is observed to increase from peripheral 
to central collisions [10]. This clearly indicates that the mechanism of particle 
production could be different in different pseudorapidity regions. In the year 
2005, a unique opportunity to investigate the system-size dependence of global 
observables occurred when Cu-|-Cu coUisions were produced at RHIC. In hght 
of the earlier results of photon multiplicity scaling with (A^part) [12,13] at for- 
ward rapidity, one could make direct comparison of the observables (A^ and 
dN^/drj) for Cu-|-Cu and Au-|-Au collisions for systems having similar values 
of (ATp^rt). 

In this paper wc present the first measurements of photon multiplicity distri- 
butions at forward rapidities in Cu+Cu collisions at ^snn = 62.4 and 200 GeV 
and Au-|-Au collisions at ^/s^^ = 200 GeV from the STAR experiment [14] 
at RHIC. The results from Au-|-Au coUisions at s/skn — 62.4 GeV were re- 
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ported in Refs. [12,13]. The photon multiphcity measurements are presented 
for various colhsion centrahty classes and are compared to corresponding re- 
sults for charged particles. The photon production is dominated by those from 
the decay of 7r°s [12]. HIJING [15] calculations indicate that about 93-96% of 
photons are from inclusive 7r° decays for the ^snn and rj range studied. 



2 Experiment and Analysis 

The STAR detector contains several subsystems which measure hadronic and 
electromagnetic observables at forward rapidity [14]. The main subsystem used 
in the present analysis is the Photon Multiplicity Detector (PMD) [17]. Pho- 
tons are detected using a highly granular preshower PMD located —5.4 m 
from the center of the Time Projection Chamber (TPC), the nominal colli- 
sion point, along the beam axis. The measurements are carried out within 
the pseudorapidity region of —3.7 to —2.3 at ^snn = 62.4 and 200 GeV for 
Au+Au and Cu+Cu collisions. A minimum bias trigger is obtained using the 
charged particle hits from the Central Trigger Barrel (CTB), an array of scin- 
tillator slats arranged in a barrel around the TPC, two Zero Degree Hadronic 
Calorimeters (ZDCs) located ±18 m from the center of the TPC, and two 
Beam-Beam Counters (BBCs) [18]. A total of 307k, 334k, 289k and 330k min- 
imum bias events for Au+Au 200, Au+Au 62.4, Cu+Cu 200 and Cu+Cu 62.4 
GeV collisions, respectively, were analyzed. These events have a collision ver- 
tex position less than 30 cm from the center of the TPC along the beam axis. 
The centrality determination in this analysis uses the uncorrected multiplicity 
of charged particles in the region \r)\ < 0.5, as measured in the TPC. The av- 
erage number of participating nucleons is obtained from Monte Carlo Clauber 
calculations [19]. The (iVpart) values corresponding to various percentages of 
the cross section for Au+Au and Cu+Cu coUisions at ^snn = 62.4 and 200 
GeV are given in the Table 1. 

The PMD consists of two planes (charged particle veto and preshower) of ar- 
rays of cellular gas proportional counters. A lead plate (3 radiation lengths 
thickness) is placed between the two planes and serves as a photon converter. 
The sensitive medium is a gas mixture of Ar and CO2 in the ratio 70:30 by 
weight. There are 41,472 hexagonal cells in each plane, which are located in- 
side 12 high voltage insulated and gas-tight chambers called supermodules. 
To each supermodule is applied a negative voltage of 1400 V as the operat- 
ing voltage. A photon traversing the converter produces an electromagnetic 
shower in the preshower plane, leading to a large signal, spread over several 
cells. In contrast, a charged particle's signal is essentially confined to a single 
cell. The photon conversion efficiency studied from simulations is found to 
increase with increasing photon energy (E^) up to 1 GeV and then saturate 
for higher energies. The typical values for the three radiation length converter 
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Table 1 

Average number of participating nucleons ((A'part)) for various collision centralities 
for Au+Au and Cu+Cu collisions at ^snn = 62.4 and 200 GeV. 
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are observed to be around 70% for Ej = 0.2 GeV and 95% for = 1 GeV. In 
the present analysis, only the data from the preshower plane are used. Further 
details of the design and characteristics of the PMD are found in Ref. [17]. 

The analysis of the data from the PMD involves the following: (a) event se- 
lection, (b) ccll-to-ccll gain calibration, and (c) reconstruction or extraction 
of the photon multiplicity. The cell-to-cell gain calibration is done using the 
ADC distributions of isolated cells (cells with six neighbouring cells having 
zero ADC). The ADC distribution of an isolated cell is treated as the re- 
sponse of the cell to charged particles, corresponding to a minimum ionizing 
particle (MIP) [17]. For most of the cells this response follows a Landau dis- 
tribution. We use the mean of the ADC distribution of isolated cells to esti- 
mate and correct the relative gains of all cells within each supermodule. The 
cell-to-cell gain variation is between 10-25% for different supermodules. The 
extraction of the photon multiplicity proceeds in two steps involving clustering 
of hits and photon-hadron discrimination. Hit clusters consist of contiguous 
cell signals. Photons are separated from charged particles using the following 
conditions: (a) the number of cells in a cluster is > 1, and (b) the cluster 
signal is larger than 3 times the average MIP response of all isolated cells in 
a supermodule. The choice of the conditions is based on results of detailed 
simulations [12,13,17]. The number of selected clusters, called '-f—like clus- 
ters (iV-j,_iike), in different supermodules for the same rj coverage are used to 
evaluate the effect of possible non-uniformity in the response of the detector. 
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Fig. 1. (color online) Top panel: Photon reconstruction efficiency (e^) (solid sym- 
bols) and purity of photon sample (/p) (open symbols) for PMD as a function of 
pseudorapidity (r/) for minimum bias Au+Au and Cu+Cu at ^/snn = 200 GeV. 
Bottom panel: Comparison between estimated e-y and /p for PMD as a function of rj 
for minimum bias Au+Au at ^snn = 62.4 GeV using HIJING and AMPT models. 
The error bars on the AMPT data are statistical and those for HIJING are within 
the symbol size. 
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To estimate the number of photons (iV^) from the detected iV^,-iikc clus- 
ters, we use the photon reconstruction efficiency (e^) and purity (/p) of the 
7-hke sample defined [12,20] as = N]if/N^, and /p = N]if/N^ -like, re- 
spectively. A'^^ ' is the number of photon clusters above the photon-hadron 
discriminator threshold. The photon multiplicity in data is then obtained as 
— (/p/e-y)A^-y_iike [12,17,20], where the ratio /p/e-y is estimated from simu- 
lations as described below. Both and fp are obtained from a detailed Monte 
Carlo simulation using HIJING version 1.382 [15] with default parameter set- 
tings and the detector simulation package geant [21], which incorporates the 
full STAR detector framework. In our previous work [13] it has been shown 
that HIJING reproduces the Nch/Nj ratio in Au-|-Au collisions at ^snn = 62.4 
GeV. For estimation of in simulations it may be important to know the in- 
clusive photon distribution. Due to lack of experimental measurement of 
the inclusive photon px distribution at forward rapidity it assumed that they 
are similar to those from HIJING model. In order to investigate the possible 
differences, e-y and fp are also obtained from a detailed Monte Carlo simu- 
lation using AMPT model [16] with default parameter settings. The AMPT 
model is a multiphase transport model which includes both initial partonic 
and final hadronic interactions. The differences between the and fp values 
estimated using the two models are less than 5%. This difference is attributed 
to systematic errors on A'^. Both and fp can vary with pseudorapidity and 
centrality. The e-y and fp for minimum bias Au-|-Au and Cu+Cu at 200 GeV 
are shown in top panel of Fig. 1. The photon reconstruction efficiency (which 
includes the detector acceptance corrections) varies from 30% at ?7 = —2.3 to 
60% at 77 = —3.7 for all collision centralities obtained from simulations for 
Au+Au and Cu+Cu coUisions at 62.4 [12,13] and 200 GeV. The purity of 
the photon sample is nearly constant as a function of f] and varies between 
40% and 60% for Au+Au and Cu+Cu colhsions at 62.4 [12,13] and 200 GeV. 
Both and fp show slight variation with system-size. The rj dependence of 
the e-y reflects mainly the varying detector acceptance between rj — —2.0 and 
77 = —3.0. There is also a small effect on the r) on e-y due to varying particle 
density as a function of 77. This is effect is already reflected in the compar- 
ison of e„, values between Au+Au and Cu+Cu. The fp values by definition 
are not affected by detector acceptance. The bottom panel of Fig. 1 shows a 
typical comparison of estimated e-y and fp using HIJING and AMPT models 
for Au+Au minimum bias collisions at 62.4 GeV. The differences are within 
5% level. The systematic errors arc discussed below. The lower limit of photon 
Pt acceptance in the PMD is estimated from detector simulations to be 20 
MeV/c. 

The systematic errors for photon multiplicity (A^-y) are due to [12,13] (a) un- 
certainty in estimates of e,y and fp values arising from splitting of clusters, 
the choice of photon-hadron discriminator threshold and choice of different 
event generators for their estimation and (b) uncertainty in iV-y arising from 
the non-uniformity of the detector response (primarily due to cell-to-cell gain 
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Fig. 2. (color online) Event-by-event photon multiplicity distributions (solid circles) 
for Au-|-Au and Cu-|-Cu at ^snn = 62.4 and 200 GeV. The distributions for top 
0-5% central Au-|-Au collisions and top 0-10% central Cu-|-Cu collisions are also 
shown (open circles) . The photon multiplicity distributions for central collisions are 
observed to be Gaussian (solid line). Only statistical errors are shown. 



variation). The error in due to (a) is estimated from Monte Carlo sim- 
ulations to be < 16% for all systems and beam energies studied. It is fairly 
independent of collision centrality. The error on due to (b) is estimated 
using average gains for normalization and by studying the azimuthal depen- 
dence of the photon density, in an rj window in the detector, and is found to 
be ~ 10%. The total systematic error in A^^ is ~19% for both central and 
peripheral collisions and is similar for Au+Au and Cu+Cu at 62.4 and 200 
GeV. The statistical errors are small and within the symbol size for the results 
shown in the figures. 
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Table 2 

Gaussian fit parameters for photon multiplicity distributions for —3.7 < rj < —2.3 
for central Au+Au (0-5%) and Cu+Cu (0-10%) at = 62.4 and 200 GeV. 



Collision Type 


(A^part) 




0"^ 


Au+Au 62.4 GeV 


347.3 


252 


30 


Au+Au 200 GeV 


352.4 


582 


52 


Cu+Cu 62.4 GeV 


96.4 


73 


13 


Cu+Cu 200 GeV 


99.0 


140 


26 



3 Multiplicity Distributions 



Figure 2 shows the photon multiplicity distributions for minimum bias Au+Au 
and Cu+Cu collisions at 62.4 and 200 GeV. The distributions for both ener- 
gies and colliding ion species show a characteristic shape with a rise at small 
multiplicity owing to peripheral events. This indicates the probability of occur- 
rence of peripheral collisions is higher. This rise is followed by a near plateau 
region with increasing photon multiplicity. This region is more prominent for 
Au+Au than Cu+Cu. It corresponds to mid-central collisions. There is a fall- 
off region in the distributions for the most central coUisions. The shape of 
the distribution in the fall-off region is governed by intrinsic fluctuations in 
the measured quantity and on the limited acceptance. Also shown in Fig. 2 
are event-by-event photon multiplicity distributions for central Au+Au (0- 
5%) and Cu+Cu (0-10%) at = 62.4 and 200 GeV. The sohd hues are 

Gaussian fits to the data. The fit parameters are given in Table 2. 



Figure 3 shows the pscudorapidity distributions of photons measured in the 
PMD for various collision ccntralities in Au+Au and Cu+Cu at ^Snn = 62.4 
and 200 GeV. As expected, the photon yield increases with decreasing The 
photon multiplicity is found to increase from peripheral to central collisions. 
Comparisons to HIJING calculations for central (0-5% for Au+Au and 0- 
10% for Cu+Cu) and 30-40% mid-central collisions are also shown in the 
figure (solid curves). The HIJING results are in reasonable agreement with 
the data for both beam energies and colliding ion species. Similar conclusions 
are drawn for other central! ty classes as well. 
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Fig. 3. (color online) Photon pseudorapidity distributions for Au+Au and Cu+Cu 
at y^SNN = 62.4 and 200 GeV. The results for several centrality classes are shown. 
The solid curves are results of HIJING simulations for central (0-5% for Au+Au 
and 0-10% for Cu+Cu) and 30-40% mid-central collisions. The errors shown are 
systematic, statistical errors are negligible in comparison. 



4 Scaling of Photon Production 



4-1 Scaling with (A^part) 



The scahng of particle multiplicity with (A^part) indicates the dominance of soft 
processes in particle production at RHIC, whereas scaling with average num- 
ber of binary collisions ((A^bin)) indicates the onset of hard processes (pQCD 
jets). The PHENIX experiment first showed that at mid-rapidity, the charged 
particle production scales with a combination of (A^part) and (A^bin) [4], indi- 
cating significant contribution of hard processes in particle production. The 
PHOBOS experiment showed that such scaling has a pseudorapidity depen- 
dence [10]. At mid-rapidity {\r]\ < 1) particle production scales with a combi- 
nation of (A^part) and (A^bin); for the range 3 < |?7| < 3.4 it scales with (A^part); 
and for the region 5 < I?;! < 5.4, the particle production per average number 
of participating nucleon pair decreases with increasing (A^part)- 
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Fig. 4. (color online) Top panel: The number of photons divided by (A^part)/2 as 
a function of average number of participating nucleons for Au+Au and Cu+Cu at 
^snn = 62.4 and 200 GeV for —3.7 < rj < —2.3 . Errors shown are systematic only 
and include those for (A'part)- Results from HIJING are shown as lines (solid for 
Au+Au and dashed for Cu+Cu). Bottom panel: Same as above, for both photons 
and charged particles from PHOBOS [10]. 
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Figure 4 (top panel) shows the variation of photon multiphcity per average 
number of participating nucleon pair with (A^part) for Au+Au and Cu+Cu at 
62.4 and 200 GeV within the range —3.7 < t] < —2.3. We observe that within 
the systematic errors, the photon multiplicity scales with (A^part) at forward 
rapidities. This indicates that the photon production at forward rapidities is 
due to soft processes. For collisions with similar (A^part), the photon multiplic- 
ity is similar for Au+Au and Cu+Cu at a given beam energy. Also shown in 
the figure are results from HIJING (solid lines for Au+Au and dashed lines 
for Cu+Cu). Considering the systematic errors shown, the HIJING results 
compare well with the data for most of the collision centralities studied. 

Figure 4 (bottom panel) shows the comparison of photon multiplicity per 
average number of participating nucleon pair vs. (A^part) and the corresponding 
data for charged particles from the PHOBOS experiment for the range —3.7 < 
T) < —2.3. Like photon production, the charged particle multiplicity at forward 
rapidities is found to scale with (A^part)- For similar (A^part)i the charged particle 
production in the region —3.7 < rj < —2.3 is also found to be similar for 
Au+Au and Cu+Cu at a given beam energy. The photon production per 
average number of participating nucleon pair is slightly lower compared to 
that for charged particles. A constant straight line combined fit to the charged 
particle results for Au+Au and Cu+Cu in Fig. 4 at ^snn = 200 GeV gives 3.8 
± 0.2, while that for photons yields 3.2 ± 0.1. For ^Snn = 62.4 GeV the values 
are 2.2 ±0.1 and 1.6 ± 0.05 for charged particles and photons, respectively. 
The difference may be due to the contribution of protons to charged particles 
at forward rapidity. The measurements (—3.7 < rj < —2.3) are carried out 
close to the fragmentation region, where protons play an increasingly larger 
role [13,22]. The ratio of the number of charged particles to photons in the 
range —3.7 < rj < —2.3 is found to be 1.4 ± 0.1 and 1.2 ± 0.1 for ^/s^^ — 
62 A GeV and 200 GeV, respectively. 



4-2 Longitudinal Scaling 

Previously it was reported that both charged particle [10,11] and photon pseu- 
dorapidity density [12,13], normalized by the average number of participating 
nucleon pairs as a function of 1] — ybcam, where ybeam is the beam rapidity, is 
independent of beam energy. Further, it was observed that such longitudinal 
scaling was centrality dependent for charged particles, but was centrality inde- 
pendent for photons [12,13]. Figure 5 shows the photon pseudorapidity density 
normalized by the average number of participating nucleon pairs as a function 
of ?7 — ybcami ioT Selected centralities (for the sake of clarity) for Au+Au and 
Cu+Cu at = 62.4 and 200 GeV. The ?/bcam values for 62.4 and 200 GeV 

are —4.19 and —5.36, respectively. The Cu+Cu results are shifted by 0.1 units 
in T] for sake of clarity. The solid line is a second order polynomial of the form 
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Fig. 5. (color online) Photon pseudorapidity distributions normalized by the average 
number of participating nucleon pairs for different collision centralities are plotted 
as a function of pseudorapidity shifted by the beam rapidity (—5.36 for 200 GeV 
and —4.19 for 62.4 GeV) for Au+Au and Cu+Cu collisions at ^snn = 62.4 and 
200 GeV. Errors are systematic only, statistical errors are negligible in comparison. 
For clarity of presentation, results for only four centralities are shown. The Cu+Cu 
data are shifted by 0.1 unit in 77 — ybeam- The solid line is a second order polynomial 
fit to the data (see text for details). 
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0.54 + 0.22{r] - t/bcam) + 0.23{r] - ^/bcam)^ fitted to all the data of Fig. 5. 
A fit to the ratio of data to this function for the results in the upper panel 
yields a value of 0.96 ± 0.01 and those on the lower panel yields 1.03 ± 0.01. 
The results demonstrate that the longitudinal scaling for produced photons 
is independent of colliding ion species. In addition we re-confirm that such 
scaling for photons is independent of beam energy and collision centrality as 
reported earlier [12,13]. 



5 Summciry 



Photon multiplicity distributions are measured at forward rapidity (—3.7 < 
T) < -2.3) for Au+Au and Cu+Cu collisions at ./s^ = 62.4 and 200 GeV 
using the photon multiphcity detector in the STAR experiment at RHIC. 
As expected, the photon yield increases with decreasing \r)\ (towards mid- 
rapidity), and is larger for collisions at higher energies. The photon multiplic- 
ity per participating nucleon pair is observed to be independent of collision 
centrality indicating that photon production is dominated by soft processes. 
A similar observation is made for charged particles, although their production 
is slightly higher. This slightly higher production of charged particles than 
photons could be due to the contribution of baryons to the charged particles, 
which may come from baryon transport and contribution from beam protons. 
On the other hand, photons are mainly from the decay of mesons (tt^). For 
collisions with similar average number of participating nucleons, the photon 
yields are similar for Au+Au and Cu+Cu for a given colliding beam energy. 
The photon production per unit rapidity per average number of participating 
nucleon pair vs. f] — yheam shows longitudinal scaling which is independent of 
beam energy, collision centrality and colliding ion species. 
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